Entanglement fidelity of quantum memories 
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We introduce a figure of merit for a quantum memory which measures the preservation of en- 
tanglement between a qubit stored in and retrieved from the memory and an auxiliary qubit. We 
consider a general quantum memory system consisting of a medium of two level absorbers, with the 
qubit to be stored encoded in a single photon. We derive an analytic expression for our figure of 
merit taking into account Gaussian fluctuations in the Hamiltonian parameters, which, for example, 
model inhomogeneous broadening and storage time dephasing. Finally we specialize to the case of 
an atomic quantum memory where fluctuations arise predominantly from Doppler broadening and 
motional dephasing. 
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The ability to store flying qubits in a quantum mem- 
ory (QM) is a fundamental component of many quan- 
tum communication schemes P, Numerous possi- 
ble methods for storing and retrieving qubits encoded 
in light pulses have been proposed 0, [j, [11 , and some of 
these proposals have recently been experimentally real- 
ized, achieving e.g. storage and retrieval of a single pho- 
ton on demand [a, f-nd entanglement between light 
and matter 0, Many promising candidate systems 
for QMs such as atomic ensembles [9| , arrays of quantum 
dots [l^ or NV centers in diamond [ll| can often effec- 
tively be described as ensembles of TV two-level absorbers 
coupling to the incoming qubit. We consider two inde- 
pendent ensembles each storing one of the logical qubit 
states. The absorbers consist of two meta-stable internal 
states I g) and | e) as shown in Fig.[lja), and a transition 
I 5) — > I e) is effected by the incoming photon in logical 
state q via coupling ilq . The states | g) and | e) are usu- 
ally not directly connected optically, with this transition 
often being achieved via an intermediate state | int) and 
additional control fields. For most of this paper details 
of such additional structure in the absorbing medium are 
not considered, and we assume that its effects on the 
properties of the absorbers can be subsumed into stochas- 
tic fluctuations of the coupling parameter fi^. After a 
storage time tg another control field is used to retrieve 
the photonic qubit. Dephasing may take place in the 
memory during the storage time, which usually leads to 
different couplings when writing and reading the qubit. 

Using these general assumptions and the notion of en- 
tanglement fidelity we derive a figure of merit T that 
measures how well a QM setup can preserve entangle- 
ment between a qubit undergoing the memory process 
(the memory qubit) and an auxiliary qubit. Our figure 
of merit J- is different from commonly-used quality mea- 
sures such as average fidelity Fa for a pre-defined set 
of input qubit states [l^l- This captures the ability of 
a memory to recreate the initial state of the qubit, and 
is equal to 1 if and only if the memory stores and re- 
trieves every state perfectly. However, depending on the 
application of the QM, one might not necessarily be con- 
cerned with exactly preserving the quantum state of the 



FIG. 1: (a) General level structure of an absorber in a QM. A 
photon with annihilation operator hq is incident on a medium 
of N absorbers, and excites one absorber into state | e) via an 
intermediate state | int) . (b) Schematic experimental setup. 
We consider a photonic qubit entangled with an auxiliary 
qubit produced by an EPR source. The photonic qubit is 
stored in the memory, and the amount of entanglement that 
remains after storage is measured. 



qubit. The preservation of entanglement might be more 
desirable in some quantum information processing and 
quantum communication schemes (T^ . ITq . for example 
in a quantum repeate r [l| , Q or in the cascaded gener- 
ation of graph states The entanglement fidelity J- 
also directly relates to the degree of violation of a Bell 
inequality by an EPR pair of photons, where one photon 
is stored and subsequently retrieved from the QM while 
the auxiliary qubit is directly detected as schematically 
shown in Fig. [TJb). The setup shown in Fig. [T]^b) could 
thus be used to measure our figure of merit. 

In the system outlined above the memory qubit is en- 
coded in a subspace of the overall photon Hilbert space 
TLa- The states of the memory and auxiliary qubits are 
denoted by | 0) and | 1). The Hilbert spaces of the auxil- 
iary qubit and the medium arc T-Lb and Tic respectively. 
The system has initial state po = | 0o) {4>Q \® Pc, where 
I 4>q) G T-La ®'Hb and pc is the initial density operator 
of the medium. We assume that the absorbers are not 
correlated initially and characterize a QM as a quantum 
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operation Am that acts on the photon as follows 



Am®I: I '^o) (00 H trc £(po) 



(1) 



I vac) represents the vacuum state and g = 0, 1 denotes 
the logical state of the qubit (the underline distinguishes 
states in Ha from memory qubit states). This annihila- 
tion operator can be written as 



where £ is a LiouviUian operating on states in Ha <8i He, 
and I is the identity operator on states in Hb- We note 
that in the work of A. K. Ekert et al. a quantum chan- 
nel for qubits A is characterized by considering the action 
of the operator A (Xi I on two qubit states. The superop- 
erator Am <X) I preserves entanglement for all two-qubit 
states only if Am is unitary (the converse is well known 
[r§|V This can be seen by using a Kraus decomposition of 
Am- We find that for at least one initial two-qubit state 
I 0e) the application of a non-unitary Am will result in a 
mixed state. Purification of {Am ^ I)(| 0e) {(f'e |) results 
in the introduction of an extra ancillary system, with 
which the memory qubit is entangled. By monogamy of 
entanglement [l^, [lo, HH , the entanglement between the 
memory and auxiliary qubits decreases. 

Motivated by these observations we write a QM entan- 
glement fidelity as follows 

^(Am) = min |(0o I ul [(Am ® I)(| M (0o I)] Um \<Po)} - 

(2) 

The quantity inside the braces is the entanglement fi- 
delity for the process Am applied to the state 
tr_B[| 0o) (00 |] (trs denotes the partial trace over Hb)- 
The entanglement fidelity was introduced as a measure to 
characterize how well entanglement is preserved by such 
a process in [l^] , and detailed discussions of its properties 
can be found in [l^, [H, . Since the standard defini- 
tion of entanglement fidelity (2^ 1 measures preservation 
of state as well as entanglement we include a unitary U m, 
which acts on Ha, to allow for evolution of the photon 
that would not decrease the entanglement present. This 
unitary is chosen to maximize !F, and thus describes an 
optimized storage process to which Am is compared in 
the same way that gate fidelity measures the success 
of a quantum gate. We also minimize over all pure two- 
qubit input states so that is a property only of the QM 
that uses the worst-case scenario as a measure of its suc- 
cess. The QM Am (and hence the LiouviUian C) consists 
of a read-in process, a period of storage, and a read-out 
process that retrieves the photon on demand a time ts 
after read-in. Note that more sophisticated choices for 
Um conditional on the outcome of measurements on the 
state of the QM after retrieving the photon might enable 
further improvement of J^. However, such schemes are 
difficult to realize experimentally and are not considered 
in this paper. Thus if .7-" = 1 wc have that Am preserves 
entanglement between the qubits, but the final and initial 
states of the photon may be deterministically different. 
The representation of the QM with Am illustrates that 
for < 1 the memory process will not be unitary. 

We now consider the photon and its interaction with 
the ensemble of absorbers. We define the annihilation 



dk5,(k)ak^> 



(3) 



where ak.A, destroys a photon with polarization Xq and 
wavevector k. The mode functions <?<j(k) are normal- 



ized. 



Oqq' 



and for simplicity we have assumed 



that each logical state has an associated single polariza- 
tion Xq . The absorbers are initially in the collective state 
I C) = I gi, . . . , (?Af); and are assumed to coherently couple 
to the photon during the whole of the read-in and read- 
out processes. The Hamiltonian for the read-in interac- 
tion between the photon in state | g) and the j"^ absorber 



is given by H, 



{a 



a^g + H.c), where a, 



U) 
eg 



I e)j- {g \ . During storage each absorber evolves according 

to the Hamiltonian iJg'^ = Sq^j{t)cf'^e , with Sq^j(t) some 
time-dependent detuning. The read-out interaction of 
the photon in logical state | g) with the absorber is mod- 



eled by the Hamiltonian iJ, 



q,j^qjaii^ +U.C.), with 
couplings ^Iqj- The dependence of the operators aqj and 
bqj on the absorber reflects the fact that due to motion 
each absorber will in general couple to a slightly differ- 
ent mode. We assume that an appropriate choice of con- 
trol held can restrict this effect to a phase 5^°-, so that 
flgj- ~ ciq cxp (i(5^°- ), and similarly for the output photon 

mode bqj — bq exp (iSq^'j)- The read- in and read-out pro- 
cesses are assumed to require a time tp each. In general 
the couplings ilqj and ^Iq j will depend on time t. In the 
following we assume a simple time dependence where the 
magnitude of the read-in (read-out) coupling is switched 
on to a constant value for the time tp that maximizes 
storage (retrieval), then switched off. For simplicity we 
also let \^q,j\ = |rigj|Vj - the generalization to differ- 
ent couplings is straightforward. Inhomogeneous broad- 
ening can furthermore lead to phases linearly increasing 
with time, and during storage some additional dephas- 
ing can occur. As a result of these assumptions we write 



nqj = Kgje'^-^'-^*) and n 



dynamics. The parameters Kqj, Kqj, Mqj, S\ 



from the 
, S':!'] and 



operator for the photon in state | g) 



a' vac > 



fq.j(ts) appears as a result of eliminating Hg ^ 

fq.j(ts) are all assumed to be real normally-distributed 
stochastic variables with respect to the storage medium. 
For instance Kqj is broadened around a mean value Kq 
by a width 'WK,q and so on. 

To obtain an analytical expression for J^, we first note 
that if Kqj = KqVj the system reduces to a two-level 
problem, and the evolution during read-in can be solved 
exactly. To this end we rewrite the read-in Hamiltonian 

Hj,'^ = ^qAe'"^' + e'^'*(e<"'* - l)]a,^,<7g) + H.c. 

where and treat the term containing the fluctuation J^^"* in Kq 



perturbatively up to second order, and similarly for Hq-'^ . 
Since any mean broadening could be corrected for, we 
assume that Kg = Alq ~ for simplicity. The general 
initial normalized photon and auxiliary qubit state can 
be written as | (/)o) = a | QO) +/? | 01) +7 1 10) +77 1 11). For 
each component of | (j)o) the evolution operator U accord- 
ing to J2j Hq^^ and Hq^^ can be used to calculate the 
final wavcfunction of the system at time t f = ts + 2tp. 
Averaging over the ensemble similarly to |25j allows us 
to rewrite Eq. ^ as 



Path 1 



Am «) I : I 0o) (00 I 



?7(|0o) (</>o|® |G) (GDiJi 



(4) 

where ((...)) denotes averaging over the stochastic 
Hamiltonian variables then tracing out the memory. 
Since | ^o) is normalized T can be calculated by minimiz- 
ing over a single parameter X = jaj + in Eq. 
This results in 

T = {A2((|6o|'))+2Ao(l-Xo)Re{((6o6I))} 

+ (1-Ao)2((|&i0)}, (5) 

where Xq is the value of X that achieves the minimization 
in Eq. ([2]), and bq is the amplitude of the final output 
photon in logical state q. Differentiating !F with respect 
to X gives a minimum of 



((hO)-Rc{((M0)} 

((i6on)-2Rc{((Mt))}+((i^in)^ 



(6) 



but if this value lies outside [0,1] then X^ = or 1. 
Applying second-order perturbation theory to the read- 
in and read-out processes as previously described gives 
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M,q) 



m{Ri + wiq) 
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N 



iiL^l)_e-K.,+<,+-/,,(t.)') 



Nil+Wlq)^ 



(7) 



Ro{{{bobi))}= n 

9=0,1 



-K,,+«'b.,+«'/,,(*=)')/2 



1 



9=0,1 



(8) 



where zDk,^ = w^^q/Kq, and 6 = (1 -|- ^ -I- iw^^q)/ (1 + 
""^K,!?)^- We see that T decreases exponentially in Wa.q, 
Wb,q and Wf^qits), and also decreases as both ^ and 
^xq/^q increase {x = K,M). Due to the factors of 
1/A'^ appearing in these latter terms, it is the exponential 
terms that will dominate for large TV. Let us also note 
that to obtain maximum absorption and emission we set 
tp — 7r/2(K^7V)i/2, so the terms containing w^, ,j could 
alternatively be seen to depend quadratically on tp. Fi- 
nally, we observe that sufficient conditions for J- 1 arc 
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FIG. 2; Experimental method of measuring T requiring stor- 
age of one logical state only. A source S produces a separable 
pair of photons, so that photon 1 is stored in the QM and pho- 
ton 2 enters the pulse shaper (PS). The photons interfere at a 
beam splitter (the PS includes a time delay), and coincidence 
measurements at detectors Di and D2 are made. 



that Wa,q,Wb,q,Wf^q{ts) < 1 aild U)K,q,UJM,q < K^, with 

the latter becoming less important as ^ oo. 

The value of Xq represents the class of states that 
achieve the minimum required in Eq. ([2]). To illustrate 
this let us consider some special cases, (i) If the states | 0) 
and 1 1_) of the photon are absorbed and emitted in the 
same way (bo = bi), then evaluating Xq gives an indeter- 
minate answer, reflecting the fact that J- is minimized by 
several choices of | (j^o). Evaluation of J- in this situation 
gives a value = (( |6o| )). (ii) If state 1 1_) is perfectly 
stored, but state | 0) is not stored at all, then bg = 0, 
Xq = 1, and ~ 0. We also compare our measure with 
the previously-defined fidelity Fa ■ If entanglement is pre- 
served i.e. = 1 then Fa = 1 if and only if the output 
photon has the same mode function as the input photon. 
In the case where the photon is stored and emitted with 
100% probability, but becomes completely decorrelated 
with the auxiliary qubit Fa could vary between and 
1 depending on the spatial mode function of the output 
photon, but T ~ 1/2. 

We now describe an experimental setup (shown in 
Fig. [2]) that, assuming case (i) above holds, would 
allow us to measure !F. After read-out but be- 
fore the beam splitter (BS) the state of the photons 



will be Pi, 



it 



E„JPrn(6r)Mvac) (vaclfori 



Po I vac) (vac I japs, where | vac) denotes the vacuum, 

aps is the annihilation operator for the mode of photon 
2 after the pulse shaper (PS), and |^™*| with m > 1 

is the set of annihilation operators corresponding to the 
eigenmodes of the state of photon 1. The eigenvalues are 
in descending order pi > P2 > ■ • • and po is the prob- 
ability of not retrieving the photon on demand. Noting 
that most detectors cannot resolve photon number, the 
probability of obtaining a click in one of the detectors Dj 
[j = l,2)isP,- = S^^iP„(l+0„)/4+po/2, and of a de- 
tection in both Di and D2 is P12 = J2m=i Pm{i-Om)/2, 
where Om = \ {va.c\aps{bm^)^ | vac) p is the overlap of 
the field modes after the BS. Both the minimum value 
of P12 and the maximum value of Pi + P2 are obtained 
when aps = 6°"' i.e. when the mode of photon 2 is pre- 
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(a) 



(b) 
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FIG. 3: The entanglement fidelity of a Raman QM with (a) 
C = 4.49 X 10"", (b) C = 2.25 x 10"^^ In both cases the 
photon bandwidth 5p = 0.1 A, with A = lO^^s 
level splittings used are | g) — > | int) = 5 x lO'^^ 



s 



^ . The atomic 
^ and 



I int) = 3.5 X lO'^'^s \ and the ensemble consisted of = 10* 
atoms. 



ciscly the dominant mode of photon 1 and for this setting 
Po + Pi = Pi + P2 — Pi2- Hence by tuning the PS the 
dominant mode of the memory photon can be found ex- 
perimentally. This tuning then corresponds to the Um 
that maximizes T as in Eq. ([2]). We can then deduce 
Pi = (( l^ol )) by removing the beam splitter and mea- 
suring the probability of the memory photon not being 
re-emitted on demand. Therefore !F can be deduced. 

We conclude our analysis by applying the fidelity mea- 
sure J- to a, specific memory setup. We determine J- for 
a QM for one single photon state based on off-resonant 
stimulated Raman scattering in an ensemble of A-atoms 
[2^. The atoms each have mass M and temperature T, 
and have the same internal level structure as the general 
absorbers considered in Fig. [H^a) . The photon is incident 
on the ensemble and excites the | g) ^ | int) transition. 
A control field drives | int) ^ | e) and stores the photon 
as a collective excitation in the ensemble. The probe and 
control fields are assumed to co-propagate with carrier 
wavevectors of magnitude kp and kc respectively. Re- 
trieval of the photon is achieved by applying another 
control field a time tg after read-in. We assume that 
the probe and control fields are both far-detuned (detun- 
ing A) from level | int), so this state can be adiabatically 
eliminated giving a medium consisting effectively of two- 
level atoms. Therefore the main source of stochastic vari- 
ation in the coupling of the atoms to the photon arises 
from the atomic motion, which we treat semiclassically 
assuming a Boltzmann distribution for atomic velocity 



components Vj in the direction of the field propagation. 
This leads to K^j = Mgj = VjU)c/c, fqj{ts) = x^j/c 
and widths given by WK.q ~ wm.q = i-^cC^^^, Wf^q = 
XC^^^7 where c is the speed of light, ( = ksT/Mc^, 
X = {kp — kc)cts, and lUc = ckc- In this scheme fq,j(ts) 
originates from the motion of the atoms during the stor- 
age time and Kqj and Mqj arise from the Doppler- 
shifting of the field frequencies. The parameter Kq = k is 
defined by the couplings of the photon and control fields 
to the atoms, and is assumed to be a constant. 

The amplitude of the final photon state can be cal- 
culated as for the general case, which upon substitution 
into Eq. ([5]) yields the following expression for JF up to 



(9) 



We see that the two main contributions to the decrease in 
J- are the Doppler broadening terms, which are quadratic 
in the ratio lOc^/C/k, and the storage time dephasing 
terms, which depend on x^C- This observation results 
in the requirement that xVC ^ 1 ii^ order to achieve 
< 1 for ^ 1. Fig. [3] shows the entanglement fidelity 
of the Raman quantum memory for two different values 
of C, and the expected decrease in J- with increasing x is 
observed. 

In summary we have introduced a figure of merit J- 
for a general QM based on gate fidelity and derived an 
analytical expression for it. Our calculations took into 
account stochastic fluctuations in the coupling parame- 
ters whose origin might vary for different QM schemes. 
We concluded by applying our formalism to a specific 
atomic quantum memory. 
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